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Abstract 

Chronic Fatigue Syndrome is a protracted illness condition (lasting even years) appearing with 
strong flu symptoms and systemic defiances by the immune system. 

Here, by means of statistical mechanics techniques, we study the most widely accepted picture for 
its genesis, namely a persistent acute mononucleosis infection, and we show how such infection may 
drive the immune system toward an out-of-equilibrium metastable state displaying chronic activation 
of both humoral and cellular responses (a state of full inflammation without a direct "causes-effect" 
reason) . 

By exploiting a bridge with a neural scenario, we mirror killer lymphocytes Tk and B cells to neurons 
and helper lymphocytes Thi , Th^ to synapses, hence showing that the immune system may experience 
the Pavlov conditional reflex phenomenon: if the exposition to a stimulus (EBV antigens) lasts for too 
long, strong internal correlations among B,Tk,Th may develop ultimately resulting in a persistent 
activation even though the stimulus itself is removed. These outcomes are corroborated by several 
experimental flndings. 

1 Introduction 

Chronic Fatigue Syndrome (CFS) refers to a clinical condition characterized by a persistent debilitating 
fatigue, neurological problems and a combination of flu-like symptoms (e.g. headache, tender lymph 
nodes), ranging from at least 6 months up to several years |1] HI [3l HI |5l |6] . 

The estimated worldwide prevalence of CFS is 0.4% — 1% (meaning over 800000 people in the United 
States and approximately 240000 in the UK) with a striking socio-economical impact: The average annual 
total value of lost productivity in the United States is $9.1 billion [7]. Such numbers propel a continued 
research to determine the cause and potential therapies for CFS, whose diagnosis is still symptom-based 
and whose origin remains elusive. 
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The resemblance of the CFS to a chronic form of acute infectious mononucleosis (AIM) has provoked 
investigation on whether this illness (whose etiologic agent is the Epstein-Barr virus, EBV), can prompt 
a chronic immune reaction in the body. 

In this work we try to deepen this point: By bridging between a neural system and an adaptive 
immune system, we show that an associative learning phenomenon might underlie a transition from an 
AIM to a CFS state. Our model focuses on the mutual interaction and regulation between B and T 
lymphocytes: when stimulated by a viral load (e.g. EBV antigens), they get activated (AIM phase); the 
active phase is then supposed to relax to a quiescent state once the viral load ceases. Actually, we evidence 
that, during such a relaxing stage, the collective behavior of the components of the system can yield non- 
trivial phenomena: if AIM phase takes a relatively long time to recover (with respect to the timescale 
that sets the standard immune response, say two weeks), which is in turn related to the success of the 
EBV to elude immuno-surveillance, a strong correlation between the activation of B and T lymphocytes 
can be accomplished. As a result, when the viral load has vanished and B cells remain activate for their 
memory role, T cells can also maintain high concentration levels since they have "learnt" that active B 
cells are associated to infectior[^ Therefore, as T cells display strong inflammatory properties, we may 
get a state of chronically active immune response (with CFS symptoms) despite the original infection is 
no longer in course. 

We stress that our approach, by applying basic concepts of statistical mechanics to immunology, points 
out emerging possible mechanisms leading to the development of the CFS. Accordingly, we overlook 
the details of the interactions at work in order to focus on the very key mechanisms underlying the 
phenomenej^ : Interestingly, as we will show, for this learning process to be properly fulfilled, T cells 
need to bypass the helper signal from specialized lymphocytes and this has been recently evidenced 
exp eriment ally. 

The paper is structured as follows: In Sec. 2 we provide a basic background about EBV, CFS and 
associative learning; these topics will be merged in Sec. 3, where we present our model. Then, in Sec. 
4, we show our analysis and results. Finally, our conclusions and discussions are in Sec. 5, while all the 
mathematics involved is reported in the appendix. 

2 Minimal background 

In this section we provide a basic background about the main features concerning chronic fatigue syn- 
drome, Epstein-Barr virus and classical conditioning, then, in the following sections we will merge such 
concepts to get an interpretation for the emergence and establishment of the CFS. 

Before proceeding, it is worth introducing the main agents of the adaptive immune system [8] , high- 
lighting the details crucial for our framework. 

2.1 The (adaptive) immune system 

An immune response is generally triggered by the introduction into the body of an antigen^ which may 
have either exogenous or endogenous origin (e.g. toxins, bacteria, viruses, cancerogenic cells). For 
instance, the genes of viruses that have infected a host cell can encode several proteins working as 
antigens. 

B lymphocytes are the agents of the humoral (i.e. mediated by secreted antibodies) immune response. 
B-cells can produce antibodies upon their full activation, and are divided in clones (ensembles of cells all 

^In this context there is no difference between activation Monte Carlothrough division among plasma and memory 
cells [8| or Monte Carlothrough Couthino idiotipic/anti-idiotypic internal images l9l 1101 : both signals are meant to work 
analogously. 

^In a very simplified parallel, phase transition classification in statistical mechanics evidences the existence of abrupt 
macroscopic changes occurring in the system under investigation, when varying its control parameters: although the mecha- 
nisms underlying e.g. the "ice- water" transition and the precipitation in an acid-base titration are completely different, the 
global phenomenology - described in the proper specific set of observable - behaves in the same way and a lot of mathematics 
and physics can be shared in their modeling (first of all the minimum energy and maximum entropy principles). 
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producing the same antibody). Activation requires antigen recognition as well as a signal from (antigen 
stimulated) Th cells. From activated B cells, specific for a given antigen, memory cells are eventually 
formed; these are long-life cells able to respond quickly to a following exposure to the same antigen. 

Furthermore, Jerne, in the 70s, suggested that antibodies must not only detect antigens, but also 
function as individual internal images of certain antigens and are themselves detected and acted upon. 
Via this mechanism, an effective network of interacting antibodies is formed, where network interactions 
provide a "dynamical memory" of the immune system, by keeping the concentrations of antibodies at 
appropriate levels. This can be understood as follows: At a given time a virus is introduced in the body 
and starts replication. As a result, at high enough concentration, it is found by the proper B-lymphocyte 
counterpart. Let us consider, for simplicity, a virus as a string of information (i.e. V = 1001001); the 
complementarjj^B-cell producing the antibody Igl, which can be thought of as the string V* — 0110110, 
will then start a clonal expansion and will release high levels of Igl. Consequently, after a while, another 
B-cell will meet Igl and, as this antibody never (macroscopically) existed before, it will attack it by 
releasing the proper immunoglobulin Ig2 (i.e., the string V** — 1001001). The latter is actually a "copy" 
(internal image) of the original virus but with no DNA or RNA charge inside. 

T lymphocytes are the agents of cellular-mediated (i.e. not involving antibodies but directly cellular 
mechanisms such as lysis) response. As in standard literature, we focus on T- helper cells {Th) and on 
T-cytotoxic or killer cells (Tk) which, when in their quiescent state, are referred to as T-CD8-t- and 
T-CD4-I-, respectively. Quiescent T-cells can be activated upon contact with cells which have previously 
interacted with the antigen: T-killer cells can interact with the so-called Class I Major Histocompati- 
bility Complex (MHC-I) expressed by all cells, while T-helper cells can interact with the so-called Class 
II Major Histocompatibility Complex (MHC-II) expressed only by antigen presenting cells (APC, e.g. 
macrophages, dendritic cells, B-cells). Active Tk expresses killer functions destroying infected cells, while 
active Th assists other white blood cells in immunologic processes, including maturation of B cells and 
activation of cytotoxic T cells. 

Every immune system cell is equipped to synthesize and release a variety of small molecules, called 
cytokines, that travel to other cells (both immune and not-immune) and up/down- regulate their growth; 
cytokines include interferons (IFNs) and interleukins (ILs). 

Before turning to our framework, built of by means of statistical-mechanics techniques, we stress that 
several other aspects and techniques stemmed from the fields of mathematics and theoretical physics are 
becoming available to investigate the biological world, ranging from kinetic theories [TT], to associative 
neural networks [12], to cellular automata and more [T31II1]. 

2.2 Chronic Fatigue Syndrome 

The literature on CFS is very broad with hundreds of analysis carried out and a rich collection of data, 
yet the clinical implications of such findings remain uncertain and a unifying, globally accepted, picture 
of its etiology and pathophysiology is still missing ^ [3l HI O [6] . 

Current theories are looking at the possibilities of neuroendocrine dysfunction, virus geneses, en- 
vironmental toxins, genetic predisposition, or a combination of these: Several researches suggest that 
Epstein-Barr Virus (EBV), by prompting a chronic immune reaction in the body, might cause CFS. 
Indeed, the phenomenology reported is consistent with the idea that the syndrome may follow the occur- 
rence of an infection yielding a massive immune response, which, for causes not yet completely clarified, 
may persist for long time, although the underlying infection is no longer in course. In fact, a CFS state 
is usually associated to an abnormal concentration and/or functioning of B-cells, T-cells and cytokines. 
Another interesting and robust immunological fact found in patients with CFS is an unusually high 
(more than 67%) increase of activated CD8-I- cytotoxic T lymphocytes with MHC-II activation markers 
dSl [M [m [m [S] ■ This will be a key point of our speculation. 

From a symptomatology viewpoint, fatigue is a common symptom, but CFS is a multi-systemic 
disease including even post-exertional malaise, unrefreshing sleep, widespread muscle, joint pain, cognitive 

''The dichotomy of a binary alphabet in strings mirrors the one of the electromagnetic field governing chemical bonds 
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difficulties, chronic (often severe) mental and physical exhaustion, muscle weakness, hypersensitivity, 
orthostatic intolerance, digestive disturbances and more. 

2.3 The Epstein-Bar virus 

EBV is one of the most successful viruses, infecting over 90% of humans and persisting for the lifetime of 
the person in a non pathogenic wa\|^ [2D1 121j . The infection can follow different pathways, in particular, 
it can turn in AIM (in up to 25% cases [H]) or it can simply introduce the virus in the host organism in 
a non apparent way. 

The virus aims to enter B-cells and, if successful, two outcomes are possible: In the first case the EBV 
begins a viral replication cycle (so called "lytic phase", a common feature of most viral infections), which 
induces the death of the infected cell, followed by the complete release of new virus particles, which are 
going to infect other cells; in the second case a state of latency (latent phase) is established where the 
"disguised" virus multiplies and stands by inside the cell, while no extracellular phenomena are observed, 
in such a way that no tackling by the immune system is evidenced. 

During the primary infection, the latent cycle and the lytic cycle proceed in parallel and the immune 
system addresses most of its resources to the lytic cycle of viral replication; the infection can be asymp- 
tomatic, have non-specific symptoms, or be so massive to result in AIM. The acute phase can last up to 
several months and it ceases when the lytic cycle is interrupted by the immune responses or by the virus 
itself, then, the infection becomes latent and the host becomes a Healthy Carrier. 

The possible persistence of the acute phase, despite a potent immune response against it, indicates 
that the virus has evolved strategies to elude the immune system. Among the different hypothesis, one 
has received particular attention 23 : the antigen BCRFlj^can simulate the signal produced by IL-10 
cytokines (which normally prompts leukocytes specialized against small-sized threatening agents, like 
EBV's antigens) and determine a delay in the immune response. More precisely, the signal from BCRFl 
inhibits the production of real IL-10; the lack of IL-10 polarizes the cellular immune response in the 
activation of a different kind of leukocyte, specialized in fighting against bigger-sized pathogens. 

We finally report an interesting study [24] on T-cell responses, in the cases of a relatively brief (2-3 
weeks) and of a protracted (4 months) acute phase. Although expansions of antigen-specific T-cells were 
observed in both situations, the T-cells reactivity occurred to be broad (i.e. addressed to several, both 
lytic and latent, antigens) and narrowly focused (i.e. mainly addressed to a singular antigen, the Early 
BMLFl), respectiveljQ 

Summarizing, a significant presence of antigen BCRFl can determine a delay in the immune response. 
As a result, the immune activity may take a long time for the clearance of the infection; during this time 
the concentration of Tk cells remains high and polarizes against BMLFl antigen as if an internal self- 
reinforcement has occurred. 

2.4 Statistical Mechanics of Pavlov effect 

Classical conditioning, experimentally demonstrated by Pavlov ^25] , is probably the most famous form 
of associative learning. The typical procedure for inducing classical conditioning on a subject (e.g. a 
dog) involves presentation of a neutral stimulus (e.g. bell ring) along with a stimulus of some significance 
(e.g. food). The neutral stimulus can be any event that does not result in an overt behavioral response 
from the subject. If the neutral and the significant stimuli are repeatedly paired, the subject eventually 
associates the two stimuli and starts to produce a behavioral response (e.g. salivation) even to the neutral 
stimulus alone. 

^strictly speaking, EBV is also associated to serious diseases as Burkitt lymphoma, but its incidence is irrelevant with 
respect to its average behavior, which is the typical outcome from a many-body theory as statistical mechanics 

^The BCRFl antigen is a Lytic Antigen sharing 70% of the human IL-IOR, which is the membrane bound receptor for 
IL-10, see also [2T] . 

^An investigation on the link between BCRFl and BMLFl can be found in |26| 
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From a statistical-mechanics point of view, classical conditioning can result from the interplay of 
dynamic phenomena, as early investigated in [27j . More precisely, statistical mechanics usually assumes 
that the states of interacting "components" are fast variables, while coupling among them evolves on 
much larger time scales, in such a way that, according to adiabatic hypothesis, the whole process results 
in two distinct time sectors; for instance, in neural scenario, learning and retrieval correspond to the fast 
(neural) and slow (synaptic) dynamics respectively |28) . Conversely, Pavlov phenomenon emerges when 
these two timescales are not so spread and a unique, coupled temporal evolution can be considered for 
retrieval variables and learning ones. 

To fix ideas, let us introduce a basic model which shall be exploited in the following. We consider 
two (on/off )-neurons (Ji = ±.1 {i — 1,2) connected by one synapse J — ±1, so that {cTi, J} G {—1,1}. 
The characteristic time for the relaxation of the two neurons is the same and denoted by r, while the 
characteristic time for the relaxation of the synapse is 0, with t ^ O. The time-averaged mean values 
of these three components are •mi{t) and wit), for which {mi{t),'w{t)} G [— l,lQ This system shows 
the capacity of a dynamical learning in the following sense: consider the action of two external signals, 
(si,S2), each applied on a different neuron ((Ti,cr2). If the stimulation of both neurons happens for a 
short time namely comparable with the short timescale (i.e., t ^ r), once one signal is removed, the 
corresponding neuron stops its activity; conversely, if the two stimuli are presented for a sufficiently long 
time (i.e. t ~ 0), due to synaptic contribution, correlations within the system develop and, if one signal is 
turned off, its corresponding neuron remains active: We will refer to this dynamical feature as associative 
learning. 

Let us deepen in more technical details the emergence of such a phenomenon. At first, both signals 
Si = 1, S2 = 1 are applied to the related neurons [regime (1, 1)], consequently, the synapse can be enforced 
(according to Hebb's prescription [29]), that is 'w{t) grows in time. After a given time tg S [0,0] one 
signal, say S2, is removed [regime (1,0)]: if an associative learning is accomplished, we expect that under 
the action of si alone, the system is still able to stimulate even a2- As we are going to show, these words 
can be translated into a system of stochastic differential equations describing the evolution of the neural 
configuration. 

Let us now outline our mapping. Neurons are B and Tk lymphocytes responding to antigenic load by 
EBV and the synapses joining them is the helper lymphocyte tuning their activity. All the lymphocytes 
producing the same antibody (BCR) for the Bs, or the same TCR for the Ts, are grouped into clones, 
whose size can be in principle (e.g. under an insult) extremely large. It should be underlined that, while 
neurons and synapses are "single" cells, the interaction in the immune level acts at the level of clones, 
hence, while single lymphocytes within each clones are clearly single cells (showing discrete symmetries 
of activity ±1, as neuronal models), their clonal behavior can be described by a continuous variable. 
Now, due to the ergodic hypothesis underlying our statistical-mechanics approach, we can interchange 
the average over time with the average over the ensamble of all the N lymphocytes cF\^\j = 1,...A^, 
belonging to the clone considered. Namely, instead of evaluating the time integral of aiit) to obtain 
its typical behavior m^, we average over the ensemble of all the same lymphocytes (hence the B-clone 
under expansion for the B, or the killer or helper for the Ts). This has the advantage of turning the 
evolution described in terms of Markov chains into Langevin stochastic dynamics, whose integration can 
be accomplished easily through standard numerical techniques. In the following, we display only the 
evolution of the averages mi — {<ti), m2 — {(T2) and w — (J), described by the system of differential 

^It is worth noting that here, coherently with the statistical mechanics philosophy, we use the ergodic hypothesis, such 
that the time average of a generic lymphocyte cTi, namely m;, can be evaluated through the ensemble average, that is 
averaging over the population of the identical lymphocytes at a given time t. Clearly, within the former approach, such 
averages must be taken over a time range at least order of r to be meaningful, allowing in this way the fast relaxation mode 
to operate. 
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equations (see Appendix for more details on its derivation): 

mi + tanh [l3{wm2 + Si)] , (1) 

K7 + tanh (/3mim2) , (2) 

7712 tanh [/3(wr7ii + S2)] . (3) 

The randomness in the stochastic evolution is ruled by /3 G Rq , which encodes the level of noise in the 
system such that for /3 = the dynamics is completely random (coherently the observables average to 
zero as they are symmetrically distributed), while for /3 — 00 the hyperbolic tangent becomes the sign 
function and the dynamics is completely deterministic. 

Finally, we stress that the statistical mechanics model we have elaborated allows a formal picture 
of phenomena which, actually, go far beyond Pavlov's conditional reflex; more generally, it describes 
processes of associative learning which has been evidenced in different biological contexts [30] , 



dm I 

dt 
dw 

Itt 
dm 2 
dt 
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3 Reading the CFS puzzle from a Pavlov perspective 

Recalling the evidence of two time-scales characterizing the evolution of EBV primary infection (see 



and Sec. 2.3), as well as the statistical mechanics model introduced in Sec. 2.4 we want to exploit the 



concept of associative learning as a bridge between the neuronal and the immune contexts; the occurrence 



of such a learning process might be interpreted as a cause for the estabhshment of a CFS (see Sec. 2.2). 

The core of our mapping is that the typical timescale for an effector response (B/K) can be smaller 
than the AIM in the long lasting case, hence allowing the development of correlation between the two 
branches. The last point could allow them to bypass helper authorization to clonal expansion. Of course, 
the statistical mechanics approach does not provide an explanation at molecular level (being unaffected 
by molecular details), but it suggests that in this context such outflanked helpers enable the genesis of 
such correlation acting as a positive (reinforcing) point of contact among them. 

Let us summarize the evolution of the phenomenology of the EBV pathology within this perspective 
in three main phases, where we also anticipate the mapping between the antigenic load on Tk and on B 
cells respectively and the state of the two stimuli (§1,52): 

• Injection and presentation, [Regime (0, 0)]. The virus starts the Lytic Cycle and infects permissive 
cells, producing antigens. The infected cells (via MHC-I) and the APCs (via MHC-II) present the 
processed antigen for specific T lymphocyte's recognition. A resting CD4-I- T-cell must be triggered 
via MHC-II present on either an APC or a (antigen-specific) B-cell (pink arrows in Figure [T]) to 
become an activated Th- 

• Activation and response, [Regime (1, 1)]. As the virus load overcomes the low dose tolerance of the 
lymphocytes, B-cells become fully activated by Th (green arrow in Fig. [I]), so that they can now 
start to differentiate into memory cells and plasma cells producing antibodies; Th cells also provide 
the second signal, indirectly via APC, for a resting CDS-f T-cell. A resting CD8-f T-cell must be 
triggered by an APC via MHC-I (first signal, pink arrow in Fig. [T]) and by Th (second signal, green 
arrow in Fig. [T]) in order to become a Tk- 

The activated Tk can now operate directly (yellow arrow in Fig. |2| on APC that is still in contact 
with the Th, or that has been previously instructed. Meanwhile, the B-cell starts to produce 
antibodies (ebv*); such antibodies, being actually new antigens for the immune system, undergo 
an equivalent recognition process (yellow pattern in Fig. [2]), so that they stimulate their specific 
B*-cells for T-cell presentation]^ 



°This standard mechanism of anti-antibodies production is due to the fact that the antibodies produced constitute a 
large concentration of proteins seen as anomalous by the host itself. We stress that, for this mechanism to hold, we do not 
need the Jerne idiotypic cascade, but only the Coutinho internal image, which has been largely revealed experimentally 

miio]. 



6 




Figure 1: Circlets represent lymphocytes (T-cells in blue, B-cells in red); the red pentagon represents 
the Epstein-Barr Virus; irregular shapes are antigens (red), antibodies ebv* (yellow), or anti-anti-body 
ebv** (azure); the rectangular shape represents APCs or infected cells (healthy B-cells are not included). 
Different stages of the process are depicted in different colors: Infection (black arrows). Presentation 
(pink arrows). Activation (green arrows). 



• Possible Learning, [Regime (1,0)]. The immune response eventually annihilates the antigenic load, 
the viral load ceases and Tk cells, no longer stimulated, can undergo apoptosis. As for B cells, they 
maintain a certain degree of stimulation (memory) even though of different nature with respect 
to the original one. In fact, as mentioned in Sec. |2.1[ after activation, B* cells secrete antibodies 
(ebv**) that, being "complementary of the complementary" [31], resemble the original virus (ebv); 
as a result, they may act as signals themselves, that is, the antibodies ebv** sustain the stimulation 
of B cells 

Now, according to the duration of the co-stimulation [Regime (1, 1)] two alternative situations 
might happen: 

— Healthy Carrier State, HCS. If the Lytic Cycle has been interrupted within a relatively short 

more simplified description would require only the presence of B-memory cells for providing signalling to Th, skipping 
any discussion on memory generation in B-cell network I31| . 




Figure 2: Kill and Memorize (yellow arrows). Agents are the same as in Fig. 
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Figure 3: Learning. Agents are the same as in Fig. [T] 



time ts by the immune response, no associative learning between the production of Tk cells 
and the production of B-cells is accomplished. The immune system has stored memory of 
the infection via memory cells and the EBV latency has established. The patient becomes a 
Healthy Carrier displaying specific memory healthy cells for those given antigens, as well as 
infected resting B cells in the Latent Cycle. 

— Chronic Fatigue Syndrome, CFS. The prolonged exposition to the (original) viral load (exper- 
imentally found to occur in the presence of large concentration of BRCFl antigen) can lead to 
an associative learning between Tk and B cells production. In fact, although Tk are no longer 
directly stimulated by the antigen, the active state of B cells can work itself as a surrogate 
stimulus. Namely, Tk bypass the Th signal and interact directly with the MHC-II signal 
provided by B-cells as APC. This is consistent with where it is reported that the BMLFl- 
specific CD8-I- T-cell (which should only recognize the class MHC-I) gets active bypassing the 
necessary T-CD4-I- indirect signal. This scenario would lead to a chronic activation state and 
it will be further discussed in the next SectiorF"! 

To summarize, this is our proposal for the CFS etiology: a massive presence of BMLFl antigens makes 
the clearance of the infection slow so that a long co-stimulation of BMLFl-specific B and Tk cells takes 
place. During this stage a learning process occurs making Tk cells able to detect the signal directly from 
B cells, hence by-passing the direct stimulation from the antigen as well as from Th- One can think at 
this situation as if Tk assumed both killer and helper functions, that is as if it switched to a hybrid state. 

In order to corroborate our speculation, we studied the dynamical properties of our model through 
analytical arguments and numerical simulations. 

3.1 Formalization 

In our interpretation, the viral load represents the external signal; when Si — +ocp^ the viral load is 
much bigger than threshold levels implied by low-dose tolerance [?, 32 , conversely, when Si — 0, it is 
much lower. Tk and B cells belonging to clones specific for EBV antigen play as the neurons cri and 

^''As already stressed, in our model there are no suggestions for this MHC-I/MHC-II switch, as these details of the 
interaction are not even introduced. On the other hand, the unbalanced Tj^jTn load emerging from our stochastic dynamics 
can be explained through this mechanism. In this sense, statistical mechanics, evidencing key mechanisms, provides hints 
for understanding experimental findings. 

^^The choice of this limit value is for simplifying calculations, the physics behind is essentially the same of every "high 
enough" load. We stress, however, that the value of the field, as introduced into our stochastic system, is not coupled to 
the noise, such that a high value implicitly accounts for low noise (i.e. high 
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<J2] the synapse J represents the Th cell. Indeed, the Th can influence both Tk and B lymphocytes, via 
the sub-populations Th^ and Tff^, respectively. Hence, the "synapsis" should be thought of as a proper 



combination of Th^ and Tff^, which results in a long relaxation time O (see Sec. 3.2). 

Following [?,[32|, the real size of the clone can be related to the time-average value of the representative 
spin, bounded in [—1, 1], by means of an exponential law: 



where i = 1, 2, and 



Mi — exp 



W = exp 



'70 



m,{t) + 1 



w{t) + 1 



(4) 



(5) 



being rji, i = 0, 1,2 a parameter that introduces the size of the pertaining population. These positions 
reproduce the chemical kinetics equation for the concentrations. Here we reasonably assume that the 
clones considered can be bounded by the same size N — 2 x 10^^, so that r/i = log(iV) [8]. 

Now, the AIM phase corresponds to a Regime (1, 1), where both Tk and B clones are stimulated by 
EBV antigens. The AIM phase is estimated to last from two weeks up to two months, so we choose t in 
the interval r G [14; 60] days. 

Once the viral shedding has been interrupted by the immune responses, the Tk clone is no longer 
stimulated, while the B-cell clone is still managing the memory of the infection (directly or Monte 
Carlothrough its conjugate specific antigen): This corresponds to the Regime (0, 1). 
The two-steps evolution described here is a useful schematization for the analytical approach developed in 
the Appendix for the special cases /? = and /3 — oo. More generally, the system of coupled differential 
equations ([I]) can be solved numerically for any value of /3 and in the presence of continuous signals. In 
particular, while S2 is still non-null over the whole range considered, si can be chosen as exponentially 
decaying (being ts characteristic time for vanishing) , in agreement with experimental findings |8j . 
Despite the system is well described by the stochastic dynamical equations (Eqs. Q), we can improve 
the picture by including proper terms accounting for the collective behavior due to interaction with other 
lymphocytes and immune agents, that is, we mimic both quiescence induction and internal signaling to 
apoptosis by introducing further two small negative fields |ei|, |e2| ^ 1, in such a way that the effective 
fields are si = si -I- ei and S2 = S2 + e2, respectively. We underline that the statistical mechanics reason 
for these small fields, whose effect would otherwise be negligible being them infinitesimal, is only breaking 
the gauge symmetry of the model, so to allow a quiescent state in the absence of signals. 



3.2 On the mixed synapse and timescales 

In our model the helper T-ccU plays the role of the synapse. Upon activation, helper T cells differentiate 
in two major subtypes known as Th^ and Th2 '■ Beyond other functions, the former maximizes the 
proliferation of cytotoxic CD8+, the latter stimulates B-cells into proliferation; also, they both produce 
cytokines which are aimed to their own proliferation and cross-regulate each other's development and 
activity 8j. The net result is that, once the Th response begins to develop, it may get polarized in one 
of the two directions (either Type 1 or Type 2), due to auto-amplification and cross regulation [H]. 

Now, when the two sub-populations are completely balanced and small (corresponding to a quiescent 
state) none of the two prevails so that, in our equivalent model, there is no link between B and Tk 
{w — 0) and no learning can be established among them (as intuitively Q — )■ oo). 

In general, one can assume that the characteristic timescales of lymphocyte growth are the same, inde- 
pendently of the particular type, that is, Tff, Tk and B cells require the same time r to adjust their 
concentrations responding to a signal. However, in our model, the three agents considered feature differ- 
ent degrees of complexity: while B and Tk can be thought of as homogenous populations, Th displays 
inner degrees of freedom, being the combination of the two sub-populations Th^ and Th2 ■ As a conse- 
quence, as the dynamics of J tunes the unbalance between these populations (and they grow closely), its 
relaxation time O is assumed larger than r. 
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Figure 4: Schematic representation of the three-spins system considered, where the complex structure 
of the synapse (central spin) is explicitly envisaged. 



Conversely, when one of the two prevails the synapse is onset {w ^ 0), so that B and Tk can 
(indirectly) interact (still retaining a large timescale 9 compared to the single clone one r). 
Hence, as envisaged by the scheme in Fig. |4] the central spin, playing the role of the synapse, can be 
thought of as a combination of two sub-populations. An effective way to relate the states of Th^ and Tf/^ 
with the overall state of Th is given by the following "average" , where we denote with wi and W2 the 
"magnetizations" corresponding to the two sub-populations: 

W1W2 r n 

W = [Wi + W2 + \Wi \ + \W2\ + \Wi - W2\\- (6) 

Indeed, this combination is in agreement with the immunological phenomenology and consistent with the 
model we are introducing: balanced, quiescent sub-populations correspond to w = 0, a large unbalance 
in favor of any of the two sub-population means w < 0, while when they are both stimulated we have a 
reinforcement effect {w > 0). 

The characteristic time for the response of B and Tk is taken to be r G [14, 60], since the AIM phase 
is estimated to range from two weeks up to two months, while for our synapse we explore & starting 
from 8 > 90 days (as the tunable parameter is the ratio t/8, the previous choice does not modify the 
results). Furthermore, while the signal S2 is constant, the signal Si (representing the real antigenic load) 
is taken exponentially decaying, in such a way that the effective time of its offset is ts', similarly, in the 
analytical approach in the appendix the signal is active for a time t € [0,ts], while the second regime 
holds for i S [ts,oo), with t < ts < &. In any case we find that the value of ts crucially determines the 
final equilibrium state. 

3.3 The role of the latent and lytic cycles 

In this section we want to deepen a way (typical of all herpes- viruses) which may further contribute to 
lengthen the hospitalization time of a CFS patient. In fact EBV, once the infection has been established 
in the host body, may hide away from immune recognition and opportunely switch between cycles of 
latency and cycles of lytic replication, somehow mirroring a switch between a quiescent and an active 
external stimulus (antigenic load) acting on the immune system. In fact, during latency, EBV mainly 
manages minimal tasks as inhibiting apoptosis and blocking viral lytic replication, while, during the 
lytic phase, EBV syntetizes proteins from many viral genes, allowing for nucleotide biosynthesis, RNA 
processing, viral DNA replication, etc. 

As a consequence, within the framework based on "Pavlov phenomenology" we are using to explain the 
transition from an AIM to a CFS scenario, these re-activations display significantly different outcomes in 
healthy carriers and in CFS patients. In fact, as shown in Fig. [Sj for the formers, whose infection walked 
off quickly (toward an healthy carrier final state), sequential impulsive stimuli do not have particular 
consequences, while, for the latters, whose infection has been prolonged enough to allow the correlation 
via helpers, basically each time there is an impulsive reactivation, this thwarts the natural de-learning. 
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Figure 5: Upper panel: impulsive stimulus s = Si = S2 as a function of time; Lower panel: immune 
response in terms of concentration of Tk (dark line) and of Th (bright line) for an healthy carrier (smaller 
values) and a CFS patient (higher values). 




4 Results 



The stochastic system of Eqs. ([T])-(3) is solved numerically by means of standard Runge-Kutta packages 
for Matlab. For the sake of clearness, a fully analytical solution of the system is shown in the appendix, 
nonetheless, here it is worth introducing the whole set of observables we need to consider. In particular, 
we have to deal with averages and correlations: 



= (J)r, mi = 



m2 = {<T2)T,mi2 = {aia2)T, 



moi = {J(Ti)t, mQ2 = (Ja2)r, "^012 = {Jcria2)T, 

where the average (•)!-, is meant over the time and over the ensamble due to ergodic hypothesis as 
explained in Sec. 2.4 (see also the appendix). We set our initial conditions m(0) using concentrations 
expressed in cells//iL and keeping in mind that, in a healthy body, a given clone has an incidence of 1 over 
10^ cells with respect to the whole population; by using the normal values for lymphocytes concentrations 
and the relative translation in terms of magnetizations (see Eq. |4| reported in Tab. 1, we get: 

w{0) - 
TOi(O) = -0.971 
TO2(0) = -0.938 

m(0) = { mi2(0) = mi(0) • 7712(0) (7) 
TOoi(O) =0 
mo2(0) =0 
moi2(0) = 

Notice that the initial value for the correlation mi2 as the product of the two concentrations mi, m2 is a 
useful condition to initialize the evolution (implicitly assuming un-correlation), and we will use it during 
numerical integration. 

In order to recover the two cases of HSC patient and CFS patient, as reported in Sec. |2.3[ we consider 
two different situations, corresponding to a short {ts = 5 days) and to a long {ts = 100 days) AIM phase, 
respectively. The other parameters holding for both patients are: 



r = 30, e = 90, 
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Table 1: Data for initial conditions [8 . Notice that we do not specify the relative concentrations of Th^ 
and Tf/^j'^tiils we assume that their sum Th is, at rest, zero, i.e. w ^ 0, meaning a balance for the two 
sub-populations. 



Agent 


Concentration 


% 


MM 


TO,(0) 




(cells/^L) 




(specific cells//iL) 


(adimens.) 


Th 


1000.50 


46 


0.010 


-0.815 


Tk 


413.25 


19 


0.004 


-0.971 


B 


500.25 


23 


0.005 


-0.938 



si = 1000, S2 = e-(*-*=)/80, 

£1 = -1.0, £2 - -1.0. 

We also fix the level of noise as low (/? — 6.5), while later we will discuss the case of high noise. 
Patient 1: HSC Scenario. 

As shown in Fig. [6j during the AIM phase {t < tg = 5 days), both signals si and S2 are active and, as 
responses, mi, m2 grow up, meaning that Tk and B clones are proliferating; the synapse also increases 
as both its afferent inputs are growing (correlations begin). 

As the real viral load is diminishing, Tfj- concentration decreases and finally reaches a value comparable 
with the initial one. Conversely, B clones, being still stimulated, maintain high levels of concentration. 

Patient 2: CFS Scenario. 

As shown in Fig. [7j during the AIM phase {t < tg — 100 days), the concentrations of B, Tx and Tff 
increase as a result of a viral load (si,S2 > 0)- This time, the signal on lasts long enough for 
Tk and Th to reach high levels (both ^ 10^ cells//iL) and, even when the signal is switched off their 
concentrations are much larger than the one pertaining to Patient 1. 

The outcomes for the two cases are compared in Fig. [Sj in terms of concentrations of Th and Tk. In 
particular, for Patient 1, Tk displays a maximum (-- 10^ cells/^L) at very short times and then relaxes 
to approximately 10^ cells//iL, while, for Patient 2, Tk reaches a higher maximum at longer times and 
then relaxes to approximately 10'^ cells//xL. 

Finally, when the level of noise is high, we expect that the effects due to interaction get more and more 
negligible. Indeed, in this model "low" and "high" levels of noise are not referred to a specific "critical 
value" , as this model does not break ergodicity by itself; conversely, since the noise level (3 is coupled to 
the averaged energy in the system E = (Jo'icr2), their product defines the levels, i.e., either (5E > 1 , or 
j3E < 1. For instance, the case /3 — 0.8 is shown in Fig. |9j notice that independently of the duration of 
Si, both the averages mi and m2 relax to small values. 

To summarize, according to the duration of the AIM phase, which is in turn related to the success of 



the EBV strategy (Sec. 2.3), we can get two possible scenarios. If the AIM phase is rather fast, when the 
viral load has vanished, B cells are continuously activated, while the concentration of Tk cells recovers 
normal values hence we reach an equilibrium state corresponding to an healthy carrier state. 
Conversely, if AIM phase is prolonged, a strong correlation between the active B and Tk lymphocytes can 
be accomplished; as a result, when the viral load has vanished, again B cells are continuously activated 
but Tk lymphocytes can maintain high concentrations. This is the actualization of a conditional reflex 
and, from a mathematically point of view it arises from the large correlation w stored, able to sustain 
the active status of Tk. 



5 Conclusion and Discussion 

The Chronic Fatigue Syndrome (CFS) has been studied for almost thirty years in the whole biological, 
medical and psychological world, being identified with tens of medical terms; as well, a full consensus on 
its genesis, physiopathology and treatment has not been reached yet. 
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Figure 6: (Color online) Top panel: Averages w, mi TO2 and signal Si as a function of time (zoom on the 
early regime). Bottom panel: the same averages and their correlations as a function of time, as shown 
in the legends; time is measured in days. The signal on Tk vanishes at around time ts = 5. The level of 
noise is low and fixed at /3 — 6.5. 



13 




Figure 7: (Color online) Top panel: Averages mi m2 and signal si as a function of time (zoom on 
the early regime). Bottom panel: the same averages ad their correlations as a function of time, as shown 
in the legends; time is measured in days. The signal on Tk vanishes at around time ts = 100. The level 
of noise is low and fixed at /? = 6.5. 
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Figure 8: (Color on line) Evolution of the concentrations of cells (histogram, left vertical axis) and 
of Tji cells (curves, right vertical axis) for the two patients considered in the case of small noise. 




Figure 9: (Color on line) Averages w, mi m2 and signal si (see the legend) as a function of time for 
Patient 1 (top panel) and Patient 2 (bottom panel). The level of noise is low and fixed dX fi — 0.8. 
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It is just the lack of a clear-cut picture that makes theoretical models very useful tools in order to 
get information about the possible causes of this disease. The present work aims to take advantage 
of the statistical-mechanics lenses to investigate why and how the CFS may establish. Of course, the 
model cannot simulate the whole immune system, which is much too complicated, rather, it has to 
get a compromise between simplification and inclusion of most important characters, the latter chosen 
according to experimental facts. 

Our framework is inspired by the well-known conditional reflex phenomenon in neurobiology, which, 
from a statistical mechanics perspective, can be recovered in terms of thermodynamic relaxation of 
complex systems [27j. Given two agents (e.g. two neurons or B and Tk lymphocytes) and a coupling 
(e.g. the synapse or the Th lymphocytes) joining them, provided that two agents are contemporary 
stimulated by two related signals (e.g. neural stimuli or antigens), then, even though one of the signal is 
switched off, the two agents may remain both active; this realizes the so-called " associative learning" . 
More precisely, along this analogy, we find that, if the Epstein-Barr virus infection is prolonged in time, 
Tff and B cells can reach high values of concentration, moreover, the latter have time enough to produce 
the image ebv**. As a result, when the (real) viral load has vanished, B cells are still stimulated, while Tr- 
may maintain high concentration levels since they have learnt the correlation between presence-of-antigen 
and B-activation. Interestingly, this learning process also implies that Tk can get activated by-passing 
the signal from Th cells and this has been recently found experimentally. 

Therefore, our model and our analysis suggest that the CFS, meant as a chronically active immune state, 
can arise from an associative-learning phenomenon. 

As a consequence, our results clarify why some individuals are affected by this disease while some other 
are not, why some infections may drive it (we investigated the case of EBV) while some other may not: 
Within our framework, the genesis of this chronically active state is guided by the ability of a particular 
virus (e.g. EBV) to elude immuno-surveillance for a larger timescale with respect to the standard ones 
set by the other antigens. Clearly, as EBV-infected patients may solve the infection in a wide range of 
times, the ones with the infection lasting for too long are natural candidates to develop CFS. Likewise, 
other pathogens, such as Citomegalovirus, can drive a long-term infection, hence allowing correlations 
and non-trivial outcomes. 

Furthermore, our suggestion may be exploited in future experiments in order to shed light on the etiology 
of this syndrome; indeed, at least at a theoretical physics level, unlearning processes are actually possible. 

Finally, it is worth stressing the "guide role" of statistical mechanics when modeling biological systems: 
in fact this approach, constraining the system to respect thermodynamics, can provide a working picture 
possibly inspiring experimental paths. In this sense, relying on natural and minimal assumptions, we 
evidence the emergence of a subtle role by Tk (able to bypass a standard double signal activation provided 
by antigen plus helpers), and, indeed, we found its consistency with experimental data. In fact in [33], it 
has been documented that a CD8-I- T receptors, that normally recognize MHC-1 signals can exhibit dual 
specificity recognizing also an antigen in the context of the MHC-II. 

6 Appendix: Evolution toward steady states of the system 

The system, whose dynamics is investigated Monte Carlothrough the paper, can be resumed as follows: 
(Ti and (72 are the effector cells (i.e., neurons in the neurobiology counterpart), while J12 is the helper cell 
(synapse); moreover we name Si,S2 the antigens (external fields) acting respectively on <Ji,a2- Lympho- 
cytes can be either active or quiescent ((7^ = ±1) and can be grouped into clones. Each clone is made up of 
a huge amount cells expressing the same idiotipicity and it is expected to be homogeneous. Analogously, 
in statistical mechanics, we can apply the ergodic hypothesis, and move from their time-average to their 
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ensemble average, using equivalently 



rrii 



T -, N 

a..m^-J2a^/\ (8) 



N 



N 

(9) 



where i = 1,2 and the index (j) denotes the element within the clone considered. We also stress that 
both the support T of the integral and the size N of the clone are large: In the limit of T, — > oo the 
variable m S K hence allowing a continuous description in terms of Langevin equations. 
Each of the variable experiences the same structure of (external/internal) fields, namely (calling as — J12 
to preserve the symmetry) 

(ai) = (tanh (^|3ai+la^+2 + s^^) ^ a^ + bi{a^+iai+2) , (10) 

where 

at = ^ [tanh(/3 + Sj) + tanh(-/3 + s^)] , (11) 
h - ^[tanh(/3 + s,)-tanh(-/3 + s,)]. (12) 
Overall, we have the following system 

(di) ^ (tanh(/3jCT2 + si)) =tanh(/3J(cr2) +si), (13) 
(J) = (tanh(/3ai(72)) - {(Jia^) tanh(/3), (14) 
(CT2) = (tanh(/3Jcri + S2)) = tanh(/3J(cri) + si, (15) 

where in the last passage of eq.s ( 13|15 ) we assumed the mean field approximation, namely 

(/(m))^/((m)), 

as we are interested in the average behavior. 

The time scale of J is 6 ^ r, being r the time scale of the a and we pose ^ + ^ = ^, ^ + ^ = ■ For 
the sake of clearness, while we used the milder notation with to, w in the text to lighten the notation, 
here we develop the whole theory using the brackets. The averages are assumed to evolve according to 



T- 



dt 



-((Ti)+ai + &i(Ja2), (16) 



= -(J)+tanh(/3)(aia2), (17) 

dt 
d{cJ2) 

T 

dt 



-{<J2)+a2 + b2{J<Ji), (18) 



while the correlations evolve according to 



^^^^ = ^(aia2) + ^ (01(^2) +6i(J)+a2(ai)+62(J)), (19) 

= ^(^^i> + |tanh(/?)(f72) + ^(ai(J)+5i(a2>), (20) 

= ^(^^2> + |tanh(/?)(fTi) + ^(a2(J)+fe2(fTi)), (21) 

^^^^^ = ^(J^i^2) + ^tanh(/3) + ^(ai(Ja2>+5i+a2(J^i)+?>2), (22) 
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as in standard ofF-cquilibrium statistical mechanics. 

The source of (ui) and of ((72) are si and S2, respectively, while the source for (J) is the internal correlation 
{0-1(^2) ■ 

According to the mutual value of Si and S2, such a system experiences different regimes, whose dynamics 
we are going to analyze. 

6.1 Regime si = 0, S2 = 0: No signalling. 

Langevin dynamics reduces to 

= -(J)+tanh(/3)(aia2), (23) 



dt 

dt 
d{(J2) 



-(ai)+tanh(/3)(Ja2), (24) 



dt 

rf(£7lfT2) 1 



-(a2)+tanh(^)(Jc7i), (25) 



dt 



T 



-2(ai(T2) +2tanh(/3)(J) , (26) 



= -^(Jai) + itanh(^)(a2), (27) 

= -^(Ja2) + ^tanh(/?)(ai), (28) 
at T r 

'-^^ = -^{Ja,a2) + ^t.nm. (29) 

As it is immediate to see, the global dynamics spreads over four different independent sub-dynamics, 

namely (cti) <^=^ {J<72), (^2) '^=^ {J<^i), {(^1(^2) *^=^ {J)i {J<^i0'2), whose asymptotic regime is given by 

(J) = (ai(T2) = 0, (ai) = (Ja2) = 0, {a^} = (Jai) = 0, (Jaia2) = tanh(/3). 

The general solution of the problem can be obtained coupling the four different sub-dynamics. 

By the first set we get 

'^^■^J* - 1 / n I tanh(^) 

- -eiJ) + —Q-M (30) 

= -'(aia2) + '*^(J). (31) 
at r T 

whose matrix can be written as 

1/e -tanh(/3)/e 
-2tanh(^)/e 2/e 

We can diagonalize the sub-dynamics by looking for solutions as linear combinations as 

Y{t)=a{Jit))+b{ai{t)a2it)), (32) 
and we can associate to this new variable a characteristic timescale f as 

f ^ = -y ^ f ( - (J) + tanh(^)(aia2)) (33) 
2bf 

h 

Namely we get the system 

fl- ^"j +2-tanh(/3)6 = (34) 
^tanh(^)a+ (1-2^)6 = 0. (35) 



+ ~ ( ~ (^i'^2) + tanh(/3)(J)) = -a{J) - b{aia2). 
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Let us work out t(/?): 



whose roots are 



1 2 



re cosh^(/3) 



0, 



n,2 (/?) 



e 



e 



Tecosh^(;9) 



(36) 



(37) 



Now we have to solve for a,b inY — a{J) + b{aia2) ■ We can define a — 2c^ tanh(/3), b — — c(l — ^), by 
which 



Yo{t) = co(2^tanh(/3)(J)-(l-g)(aia2)) 
Yi2{t) = ci2(2^tanh(;3)(J)-(l-g)(aia2)), 



(38) 
(39) 



on which we can fix cq,ci2 as cq = r/2Ti tanh(/3), ci2 = —1/(1 — T2/e). From eq.(36l we get 

tanh(/3) = 



re /I 1 



e 



1 2 

f T 



and we can solve for Yq, Y12: 



Yo{t) = (J) + 

yi2(t) = (^1(72) 



1 

e fi 



J 2 

T2 T 



f2 



SO we get the form 



with 



Yo{t) 

Yl2{t) 



{J{t))+Ao{(Ji{t)a2it)), 
{'Ji{t)a2{t))-A,2{Jit)), 



(40) 
(41) 



(42) 
(43) 



^0 



j_ i_ 

Ti_ 

2 1_ 

T ti 



A,2 = 



J 2 

■f2 T 

\ i - i 

\ T2 e 



By the eigenvalues found in eq.s(37l we can build the eigenvectors Vi — (Vn, V12), V2 ~ (V21, V22) as 



Vn = tanh(/3)/e, F12 



^21 = tanh(/3)/e, 1/22 



2r r 2 / ' 



1 1 1 

26 " 7 ^ 2 



being A = (1/9 + 2/t)2 - by which, finally we get 



(/9)' 
(^lW'^2(t)> 



CiViie ^1 + (72^126 *2 , 

ClF2ie"^ +(72^226"^. 



(44) 
(45) 



By the second set we get 
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= -(ai)+tanh(/3)(Jr72): (46) 
r'^^^ = -(Ja2)+tanh(;3)(ai). (47) 

whose matrix can be written as 

1/t -tanh(;3)/T 
-tanh(;3)/T' 1/t' 

Again we can write a solution in the general form Y{t) = a{ai) + b{Ja2) and label f its characteristic 
timescale such that 

dt 



a(- ((Ti) +tanh(/3)(Ja2)) + ;^6( - (Jaa) + tanh(/3)(ai)) 
= -a(<Ti) + b{Ja2), 



and write the system 



(1 - -)a+ ^tanh(/3)6 = 0, (48) 

T T 

- tanh(/3)a + (1 - = 0. (49) 

T T 



Again we can find the eigenvalues 2 iP) 



fi,2(/3) 2Vr r' VV T 
and write the general solution in the form 

, ,,, Citanh(/3) _i C2tanh(/3) _^ 

((71 (t)) = ^1 ^2, (51) 

T r 

(Ja2(i)) = -^(^-VA)e-'^-^(l-yA)e-^, (52) 
A being (1/r + l/r'f - {4/tt'){1 - tanh(/3)). 

For the two other subsystems we can proceed exactly as we did so far and verify that, even though 
on different timescales, all the observables (magnetizations and correlations) converge to zero. 
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6.2 Regime si = 0, S2 = oo: One infinite signal. 

Let us start with the following conditions: si — 0,S2 = oo, then we have ai = 0, 61 = tanh(/3),a2 = 
1, 62 = and the evolution of the system can be written as 

= -(J)+tanh(/3)(aia2), (53) 

T — ; — = 
dt 

d{cJ2) 
T — ; — = 
dt 

d{(Ji<J2) 

T = 

dt 

= -^(Jai) + itanh(/3)(a2), (57) 

5^ = _i_(j^,) + ^tanh(^)(ai) + ^(J), (58) 

d{J^ = -4(Jaia2)+^(Jai) + ^tanh(/3). (59) 

at T T T 

The presence of a signal in the second channel bridges the two sets {{ui (t)), {Ja2{t))) and {{(Ji{t)a2{t)) , {J{t))) 
such that the latter, if the system has experienced the field enough time to learn correlations, will assume 
high values: a signal in the second channel may induce a response even in the first channel. 



-(ai)+tanh(^)(Ja2), (54) 
-(^72) + 1, (55) 
-2(<Tia2) + tanh(^)( J) + (<ti), (56) 



Let us divide the system again by considering the following, natural sub-dynamics 

dt 

d{(Tia2) 



qM^ = -(J)+tanh(^)(aia2), (60) 
at 



. ■ -2(criCT2) +2tanh(/3)(J) 

+ [(ai)-tanh(/3)(J(f))]. (61) 

We observe that the term [(cti) — tanh(/3)(J(t))] is a novelty with respect to the same equations in the 
regime .si = 0, .S2 = and represents the "unlearning source": if the conditional reflex (which we are 
going to introduce) is not reinforced, it will tend to vanish. 

The source of such a conditional reflex can be found by looking at the other subsystem, namely 

d{ai) 



dt 



-(ai)+tanh(/3)(Ja2), (62) 



^ = -^(Ja2) + ^tanh(/3)(a,) 

+ -[(J)-tanh(/3)(ai)]. (63) 

T 

In fact, the new term [{J) — tanh(/3)(cri)], namely the source of the conditional reflex, is a learning term 
as it couples the slow channel ((J), {(Jia2)) with the fast one ((cri), (J(T2)). 

We can now start studying the dynamics in this regime by considering the sub-dynamics of ((J), (ci), 
(ci(72), {J<J2)), whose dynamical system reads off as 

d{ai) 



dt 



-(ai)+tanh(/3)(Ja2), (64) 



= -(J)+tanh(/3)(c7ia2), (65) 

_d{aia2) 



dt 



-2(aia2) + tanh(/3)( J) + (cti), (66) 



= -^(Ja2) + ^(J) + ^tanh(^)(<7i). (67) 
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The associated matrix can be written as 



1/r 


1 



V -itanh(^) 





1/e 

itanh(/3) 





-tanh(/3)/e 

2 

T 





tanh(/3)/r \ 




We can diagonalize the dynamics and look for solutions as linear combinations like Y{t) = a{ai)+b{J) + 
c{(J\(J2) + d{Ja2), associating to this variable its characteristic timescale f ^^Y, and proceed as for the 
former regime. 

Skipping all the calculations for the sake of brevity we report only the solution 

t_ I ^ // ^ /// ^ 

{ai{t)) = cixie ^1 +02X^6 ^2 +C3.x^e ^3 + c4Xj e ^4, 

t_ I t_ II t_ /// t_ 

{J{t)) = ciX2e ^1 +02X26 ^2 +03X26 ^3 +04X2 e ^4, 

^ / ^ // ^ /// t_ 

{ai{t)a2{t)) = 01X36 ^i+C2a;3e ^2+030:36 ^3+04X36 ^4, 

(J(i)cr2(t)) = cia;4e ^1 + 020:46 ^2+03X46 ^3+04X46 ^4, 

all the x's being the component of the following eigenvectors Vi = {xi,X2,X3,xA, V2 = (a;'i , a;,, a;',, a;',), 

_ // // // // _ /// /// /// /// 
V3 — yx-^ ,X2, X3, X3 ), Vi — yx-^^ ,X2 1X3 ,X3 ): 



Xi 
X2 
X3 
X4 



-tanh(^), 
r 

2Vt e 

-tanh(^), 

T 

1 /I 



2 Vr 



1 

e 



- VA 



= -tanh(/3), 



tanh(^), 



2 Vr 



1 

e 



= — tanh(/3), 
r 

= 

2 Vr e 
e [- 

~ tanh(/3) 
2Vt 



2x9 



1 

292 



;^tanh2(/?) 



A(- 



26^ 



22 



X-, = 



Xo — 



X4 — 



-tanh(/3), 

2\t e 
e 



tanh(^) 



2tQ 292 T 
1 

2 Vr e ' 



, ^tanh^(/3)-VA(i-^) 



being A = l/r^ + l/G^ + 2(2tanh2(/3) - l)/(Te). 

The remaining sub-dynamics of (((T2), (Jcri), (Jcricr2)) is depicted by the system 



T ; 

dt 
,d{Jai) 

T ; 

dt 

d{JaiCF2) 
dt 



-(^2) + l, 

-(Jai) + (a2), 
\,{J<Ji<y2) + -{Jcri) + ^ tanh(/3). 



T" T T 

By applying the framework previously shown several times, we obtain the solutions 

t_ 

(0-2) = cie ^1 + 1, 

(JcTi) = ci— e "1 + C2e ^2 + 1, 



(JcrifT2) =ci— e ^1 +C2e ^2 +036 ^3 +r 



- tanh(/3) 



(68) 



(69) 
(70) 
(71) 



6.3 Regime si = 00, S2 = 00: Two infinite signals. 

Let us start with the following conditions: si = 00, S2 = 00, then we have ai = 02 = 1, 61 = 62 = 0) and 
the evolution of the system can be written as 

d{J) 
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dt 
d{ai) 

T ; 

dt 
d{(J2) 

T ; 

dt 
d{aia2) 

T ; 

dt 

d{Jai) 



= -(J)+tanh(/3)((TiCT2), 
= -{ai) + l, 
= -{cr2) + l, 



= -2(<TlO-2) + (ai) + ((T2), 



dt 
d{Ja2) 

dt 



= - ^(Jai) + ^ tanh(^)(c72) + -{J), 
r' 9 r 

= - A(J^2) + ^ tanh(/3)(ai) + -{J), 



d{Jaia2) 
dt 



9 

^{J(Tia2) + ^ tanh(/3) + -(J(T2) + -(Jcti), 
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whose solutions, in complete analogy with the previously introduced methodology, can be obtained as 



{J{t)a,{t)) 
{J{t)a2{t)) 
{J{t)ax{t)a2{t)) 



{J{t)) = cre-i +{C2 + C^)^^^e- 

036"- + 1, 

(C2 + caje T + C4e ^ + 1, 

t tanh(/?) , C2 C3 , t 

(5-7) ^ ^ 

t tanh(/3) , C2 C3 , * 
cie-e + ^—^(^ + ^)e-7 + C4 

__L tanh(/3) , x _i 

2cie e + ^1 ' / (C2 + C3)e - + C4 



tanh(/3) _2t , 
C4 _ " +tanh(/3), 



tanh(^) 



tanh(/?) 
tanh(/3) 



+ cye"^ +tanh(^) + (c5 + C6)e"^ +tanh(^). 



-e ^ +056 ^+tanh(/3), 
-e~^ +C6e~^ +tanh(/3), 
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